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ABSTRACT 

The energy to desorb atomic oxygen from an interstellar dust grain surface, E^^s, is an important 
controlling parameter in gas-grain models; its value impacts the temperature range over which oxygen 
resides on a dust grain. However, no prior measurement has been done of the desorption energy. We 
report the first direct measurement of Edes for atomic oxygen from dust grain analogs. The values of 
Ades are 1660 ±60 K and 1850 ±90 K for porous amorphous water ice and for a bare amorphous silicate 
film, respectively, or about twice the value previously adopted in simulations of the chemical evolution 
of a cloud. We use the new values to study oxygen chemistry as a function of depth in a molecular 
cloud. For n = 10^ cm“^ and Go=10^ (Go=l is the average local interstellar radiation field), the 
main result of the adoption of the higher oxygen binding energy is that H 2 O can form on grains at 
lower visual extinction Ay, closer to the cloud surface. A higher binding energy of O results in more 
formation of OH and H 2 O on grains, which are subsequently desorbed by far-ultraviolet radiation, 
with consequences for gas-phase chemistry. For higher values of n and Gq, the higher binding energy 
can lead to a large increase in the column of H 2 O but a decrease in the column of O 2 . 

Subject headings: astrochemistry—ISMiabundances—ISMimolecules—ISMiatoms—ISMidust, extinc¬ 
tion 


1 . INTRODUGTION 

Although oxygen is the third most abundant element in 
the universe, finding it has proven to be challenging, es¬ 
pecially in molecular form in dense clouds. Observations 
of molecular oxygen with the Submillimeter Wave As¬ 
tronomy Satellite (SWAS) toward nearby clouds yielded 
no detection, while with the Odin satellite upper limits of 
10 “^ (referred to H 2 ) were obtained for several cold dark 
clouds that were examined (Pagani et al. 2003; Larsson 
et al. 2007). Yildiz et al. (2013) put an upper limit of 6 x 
10“® in a Herschel HIFI search toward a low-mass pro¬ 
tostar. Detection with Herschel HIFI of O 2 in Orion, via 
the rotational lines at 487, 774, and 1121 GHz, gave an 
estimated O 2 abundance of 10“® (Goldsmith et al. 2011); 
such a large value is attributed to the passage of a shock 
(Chen et al. 2014). Towards p Oph A, the estimated 
O 2 abundance is 5 xl0“® (Larsson et al. 2007; Liseau 
et al. 2012). No detection at 487 and 774 GHz toward 
the Orion Bar, an far-ultraviolet (FUV) illuminated re¬ 
gion at the interface of H H region and a dense molecular 
cloud, was reported (Melnick et al. 2012). Overall these 
observations confirmed that molecular oxygen abundance 
in dense clouds is much less than predicted by gas-phase 
models, that is, 7 x 10“® for O 2 in a cold molecular 
cloud in steady state (table 9 in Woodall et al. 2007). 
Hincelin et al. (2011) looked at gas-grain oxygen chem¬ 
istry and studied the sensitivity of rate coefficients on O 2 
abundance. Motivated by the conflict between gas-phase 
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models and observations, Hollenbach et al. (2009) con¬ 
structed a gas-grain model of the chemistry of molecular 
cloud surfaces and interiors. They found that the low 
O 2 abundances could be explained by models in which O 
was largely atomic at cloud surfaces and locked primarily 
in water ice deeper into the cloud, with smaller contri¬ 
butions from CO, CO ice and refractory grain material. 
This model and others (e.g. Garrod et al. 2008; Hasegawa 
& Herbst 1993) highlight the importance of grain-surface 
reactions in understanding abundances in the interstel¬ 
lar medium. Depletion of oxygen in silicate grains is well 
constrained, as well as in CO and N-containing molecules 
(Jenkins 2009, 2014). Oxygen is present in water, but 
gas-phase water abundance is not enough to explain the 
missing oxygen. Water is also in ice-coated grains, but 
more water ice than it is observed in ice is necessary to 
account for all of the oxygen. Some water can be in 
hydrated silicates, on grains in low Ay regions before 
thick ices are formed, or on large (> 1 pm) grains; how¬ 
ever, each of these proposed oxygen reservoirs has its 
own weaknesses (Whittet 2010). It has been suggested 
that oxygen is trapped in other types of grains, such as 
oxygen-bearing carbonaceous grains, similar to organics 
in cometary particles (Whittet 2010). 

On grains, oxygen can react with hydrogen to yield OH 
and then H 2 O via hydrogen addition reactions, or it can 
diffuse and form molecular oxygen and ozone (O is al¬ 
most exclusively converted into OH/H 2 O on the surface. 
O 2 is likely formed in the gas phase and eventually con¬ 
denses on grain surfaces). From laboratory experiments, 
we know that O 2 leaves the surface of dust grain analogs 
at ^ 30 K, but O 3 at a much higher temperature, ^ 67 
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K (He et al. 2014). O 3 can be hydrogenated on grains 
and form water (Mokrane et al. 2009; Romanzin et al. 
2011; He & Vidali 2014). In recent years several experi¬ 
ments were done on dust analogs in order to explore the 
routes to water formation via hydrogenation of O, O 2 
and O 3 (see references in Vidali 2013). How successful 
these routes are in forming water or hydroxyl products 
depends in large part on how successfully a dust grain 
keeps O on its surface. Tielens & Hagen (1982) derived 
the adsorption energy of O on grains, by assuming 
that it is held by physisorption forces that depend on 
the atomic polarizability and the dielectric function of 
the solid (Bruch et al. 2007). Making an analogy with 
adsorption of gases on graphite for which an extensive lit¬ 
erature exists (Vidali et al. 1991; Bruch et al. 2007), Tie- 
lens & Hagen (1982) derived iilb=800 K. Eventually, this 
value was adopted in codes of simulations of gas-grain 
processes in the ISM (e.g. Stantcheva et al. 2002; Garrod 
et al. 2008). One case where the models of Hollenbach 
et al. (2009) predicted high O 2 abundances (~ 10“®) was 
in regions with very high FUV fields, where desorption of 
O from grains would allow gas-phase production of O 2 . 
Recognizing that observations were in conflict with this, 
they suggested the O binding energy was at least 1200 K. 
Observations of the Orion Bar by Melnick et al. (2012) 
set a stricter limit on the O 2 abundance in a high FUV 
region, leading these authors to propose that the bind¬ 
ing energy was as high as 1600 K. He et al. (2014) used 
a rate equation model to analyze the data of O 3 forma¬ 
tion on an amorphous silicate in the temperature range 
of 30-50 K and indirectly obtained the atomic oxygen 
adsorption energy, 1760 ± 230 K. 

Here we present a direct measurement of the desorption 
energy Udes of atomic oxygen from porous amorphous 
water ice and from a bare amorphous silicate. (The des¬ 
orption energy is the binding energy E^, if there is no 
activation process, as we assume here in absence of evi¬ 
dence pointing otherwise). Then, we use the photodisso¬ 
ciation region (PDR) model of Hollenbach et al. (2009) 
to look at the consequences that this revision has on the 
chemistry of the ISM by doing a simulation with the old 
and the new values of the desorption energy. We also look 
at the effect on ISM oxygen chemistry of using a revised 
value for OH desorption energy (He & Vidali 2014). 

2 . EXPERIMENTAL 

Here we give a brief description of the experimental set¬ 
up and measuring methods; details can be found in Jing 
et al. (2013) and He et al. (2014). The apparatus con¬ 
sists of a triply differentially pumped atomic/molecular 
beam line connected to an ultra-high vacuum chamber. 
Molecular oxygen is dissociated in a water-cooled Pyrex 
tubing placed in a radio frequency (RE) cavity. The O 2 
gas flow is controlled by a mass flow controller. The 
beam enters the sample chamber through a 2.5 mm col¬ 
limator. The narrow angular spread and the triple dif¬ 
ferential pumping of the beam line minimize the gas load 
to the sample chamber and the deposition of beam parti¬ 
cles on the sample holder. The ultra high vacuum cham¬ 
ber is pumped by a turbomolecular pump, a cryopump 
and an ion pump. The background pressure is about 1 
X 10“^® torr when the system is at room temperature 
and 5 x 10“^^ torr when the sample holder column is 
cooled. When the beam enters the chamber, the pres¬ 


sure is of the order of 1-2 x 10“^° torr. The detector 
is a Hiden Analytical HAL-5 quadrupole mass spectrom¬ 
eter (QMS) that is mounted on a doubly differentially 
pumped rotatable platform. Therefore, the detector can 
be swung in front of the beam to characterize the incom¬ 
ing beam and can be placed facing the sample during 
the deposition and the temperature programmed desorp¬ 
tion (TPD). The detector itself is placed in an enclosure 
that is differentially pumped. The aperture to the ion¬ 
izer of the QMS is fitted with a Teflon cone with the 
apex facing the sample. The distance between the tip 
of the cone and the sample can be adjusted by reposi¬ 
tioning the axis of the sample via an XYZ manipulator. 
The purpose of placing the cone is to maximize the flux 
coming from the sample and minimize the flux of par¬ 
ticles coming from other parts of the sample holder. In 
this way, the signal/noise is enhanced. The procedures 
to calibrate the beam flux are described in (He et al. 
2015). When the molecular oxygen beam is about 0.34 
ML minute”^ (3.4 x 10^^ cm“^minute“^), the efficiency 
of dissociation is 42 %. Thus the beam intensities of O 
and O 2 are 0.29 ML minute”^ (2.9 x 10^^ cm“^minute“^) 
and 0.20 ML minute“^ (2.0 x 10^^ cm“^minute“^), re¬ 
spectively. 

The preparation and characterization procedures of the 
I /im thick amorphous silicate film have been described 
before (Jing et al. 2013). The sample was mounted on a 
sample holder that was cooled with liquid helium. Prior 
to a set of experimental runs the sample was taken to 
400 K for cleaning. We found that this procedure yields 
the same results as by sputtering the sample with Ar 
ions (Jing et al. 2013). The sample temperature was 
measured by a calibrated silicon diode placed behind the 
sample; a thin silver foil assures a good thermal contact 
between sensor and sample. There is a thermal switch 
between the sample holder and the liquid helium stage 
that allows the heating of the sample without warming 
up other parts of the apparatus. This assured a lower 
background gas pressure during the experiments. 

During the atomic oxygen (O) exposure, because the 
dissociation rate is not unity, there is always O 2 mixed 
with O, so ozone formation via O-I-O 2 —?> O 3 is likely. 
When O desorption is seen in the TPD, it could come 
from the bare silicate surface and/or from ozone that 
was formed on the silicate as explained above. To find 
out which one is true, the silicate surface was pre-covered 
with different coverage of ozone before depositing atomic 
oxygen. If O desorbs from the silicate surface, then as 
the pre-covered ozone approaches one monolayer, the 
amount of O desorption should decrease to zero. In the 
other case, if O desorbs from an ozone patch, then as 
the pre-covered ozone amount increases from zero to one 
monolayer, the amount of O desorption should increase 
linearly. The correlation between the coverage of ozone 
and the amount of atomic oxygen direct desorption in 
the TPD can suggest whether O is desorbing from the 
silicate surface or from ozone. It will be shown later that 
the former is true. The presence of ozone on the sili¬ 
cate surface actually helps to decrease the formation of 
O 2 and O 3 via, respectively, the 0-1-0 and the O-I-O 2 
reactions, since O and O 2 do not react with O 3 . Ozone 
functions as an obstacle in the diffusion of O and O 2 on 
the surface. The procedure to prepare the sample with a 
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layer of ozone was described in He et al. (2014). Briefly, 
the silicate surface was exposed to the dissociated oxy¬ 
gen beam and ozone was formed by O-I-O 2 —t O 3 ; then 
the surface was annealed at 50 K to desorb unreacted 
02 - During the exposure stage of ozone preparation, the 
beam flux was set to a higher value than that for the O 
deposition, so that the dissociation rate was lower. The 
beam intensity of O 2 was higher than O so that almost 
all the O was reacted. 

In dense clouds, dust grains are usually covered by ice 
mantles, with water as the main component. We also 
carried out measurements on porous amorphous water 
ice. The ice substrate was prepared with background de¬ 
position through a microcapillary doser. The thickness 
was estimated from the background pressure and the du¬ 
ration of the deposition. A 500 ML (with an uncertainty 
of 30 %) porous amorphous water ice film was deposited 
at 70 K on top of a 500 ML single crystalline water ice. 
The ice film was cooled down to lower temperatures for 
O/O 2 beam irradiation and TPD. During the following 
TPD experiments, the porous ice film was maintained 
below 70 K to prevent pore collapse in the ice films (see 
the temperature ramp in Figure 1). 

2.1. Results and analysis 

Figure 1 shows TPD spectra of oxygen from a porous 
amorphous water ice film. We deposited 480 s of O/O 2 on 
porous amorphous water ice at 48 K. After the exposure, 
the surface was cooled down to about 20 K and then 
heated up at a linear ramp rate of 0.167 K s“^ to 70 
K to desorb the atoms and molecules on the surface. 
During the TPD, the desorption rate was recorded by 
the QMS facing the sample. Ozone was formed on the 
ice and came off at about 67 K. Because ozone is fragile 
against electron impact ionization, a good fraction breaks 
up into and O"^, as it can be seen from the peaks of 
mass 16 amu and 32 amu at about 67 K. In the mass 
16 amu signal, there is a peak at 57 K, which is due 
to O direct desorption. As the sample was heated, O 
diffuses and forms O 2 , and O 2 leaves the surface upon 
formation, because of the short residence time of O 2 . 
The mass 32 amu peak centered at about 50 K is due to 
O 2 formation via O diffusion. As O is gradually used up, 
the O 2 formation rate decreases, and the rate of O direct 
desorption increases. The TPD peak area of O direct 
desorption corresponds to about 3% of a layer. 

Similar experiments were performed on the surface of 
an amorphous silicate. Compared to porous amorphous 
water ice surface, the amorphous silicate surface is flat¬ 
ter. The mobilities of O and O 2 on it are higher, and the 
formation of O 3 and O 2 are also faster. To hinder the for¬ 
mation of O 2 and O 3 , and leave more atomic oxygen on 
the surface available for direct desorption, we pre-coated 
the surface with a fraction of a layer of O 3 . O 3 does 
not react with atomic oxygen or molecular oxygen; O 3 is 
used to hinder the diffusion rate of atomic and molecular 
oxygen. To find out whether atomic oxygen is desorb¬ 
ing from bare amorphous silicate surface or ozone patch, 
measurements were performed to check the amount of O 
direct desorption from silicate surface covered with 0 . 2 , 
0.4, 0.6, and 0.8 ML of ozone. Figure 2 shows the TPD 
spectra after depositing 480 s of O/O 2 on a 0.2 ML O 3 
pre-coated amorphous silicate surface. In the figure, the 
two contributions to the mass 16 amu peak, i.e., O direct 



Fig. 1.— TPD traces of mass 16, 32, and 48 amu after depositing 
480 s of O/O 2 at 48 K on porous amorphous water ice. The heating 
ramp is 0.17 K s“^. The dashed line shows a fitting of the mass 
16 amu peak using E'des(O) = 1660 K (see text). 


desorption and fragmentation, are merged into a sin¬ 
gle broad peak. To find out the peak position and peak 
area of O direct desorption, the contribution needs to 
be subtracted out. The fragmentation pattern of is 
measured by O 3 TPD experiments and not shown here. 
The raw TPD data and that after subtracting con¬ 
tributions are shown in Figure 3. In each panel of this 
figure, the first row has the raw TPD data, while in the 
bottom row the contributions of ozone fragmentation to 
mass 16 amu have been subtracted to show the amount 
of direct O desorption. After this subtraction, the O 
direct desorption peak is centered at about 64 K and is 
clearly different from the ozone desorption peak. We also 
found that as the coverage of the pre-coating ozone layer 
increases from 0.2 ML to 0.4 ML, the amount of O di¬ 
rect desorption increases. This is because the increased 
ozone amount helps to lower the O diffusion rate. How¬ 
ever, when the pre-coating ozone is more than 0.6 ML, 
there is almost no O direct desorption peak. This indi¬ 
cates that the observed O direct desorption is from bare 
silicate surface instead of ozone ice; otherwise, the O di¬ 
rect desorption amount should be positively correlated 
with the amount of ozone pre-coating. We also checked 
that the desorption energy of O atoms is not influenced 
by the presence of ozone at low coverage. In a work 
by Bennett & Kaiser (2005), the formation of a O 3 -O 
complex was deduced from IR data in an experiment of 
bombardment of an O 2 ice at lOK by 5keV electrons. In 
our experiment we find that the mass 16 amu peak shape 
and temperature, see Figure 3, do not change as a func¬ 
tion of O 3 pre-coverage (at low coverage) suggesting that 
there is no effect of a possible O 3 -O interaction on the 
desorption on O. Furthermore, there is no reaction of O 3 
with O because of an activation energy barrier (Atkinson 
et al. 2004). 

The desorption energy of atomic oxygen from differ¬ 
ent surfaces can be calculated using an Arrhenius-type 
expression: 


dN{t) 

dt 


= —vN{t) exp( 


Rdes 

kBT{t) 


) 


( 1 ) 


where N (t) is the number of molecules on the surface at 
time t; iz is a pre-exponential factor which is standardly 
taken to be 10 ^^ s“^; E^es is the desorption energy; is 
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Temperature (K) 

Fig. 2.— TPD traces of 16, 32, and 48 amu after depositing 
480 s of O/O 2 at 48 K on an amorphous silicate that has been 
pre-coated with 0.2 ML of ozone.The heating ramp is 0.5 K s~^. 
The inset shows the temperature ramp. 

the Boltzmann constant; T(t) is the surface temperature. 
In Figure 1 the dashed line shows the fitting of mass 
16 amu trace using E^es = 1660 K. In Figure 4, the 
mass 16 amu traces in the top-right panel of Figure 3 
after subtraction of O 3 contribution are plotted together 
with a fitting using E^es = 1850 K. The energy values 
are shown in Table 1. The error bars of the desorption 
peak temperature are determined by the uncertainty in 
locating the peak positions. The width and shape of 
the TPD peak are not taken into account. Therefore, 
the error bar of desorption energy does not represent 
the width of distribution; instead, it represents only the 
uncertainty in finding the peak of the TPD trace. 

TABLE 1 

Desorption energy of atomic oxygen from porous 


AMORPHOUS WATER 

ICE AND FROM AN AMORPHOUS 

SURFACE. 

SILICATE 

Surface 

Desorption Peak Temperature 

i?des(0) 

Porous water ice 

56±2 K 

1660±60 K 

Amorphous silicate 

64±3 K 

1850±90 K 


3. SIMULATIONS AND ASTROPHYSICAL 
IMPLICATIONS 

In order to explore the implications of a higher oxy¬ 
gen desorption energy, we have used the PDR model 
of Hollenbach et al. (2009) which computes the chemi¬ 
cal and thermal structure of a molecular cloud edge ex¬ 
posed to FUV radiation from either the interstellar ra¬ 
diation field or from nearby hot stars. The important 
input parameters are the gas density, nn, the strength 
of the FUV radiation field, Go (where a value Gq = 1 
corresponds to the FUV radiation field in the diffuse in¬ 
terstellar medium), the abundances of important atoms, 
and the properties of interstellar dust grains. With these 
inputs, the model self-consistently computes the gas and 
grain temperatures, and the abundances of more than 60 
species in the gas phase and in ice mantles on grain sur¬ 
faces, starting from the cloud surface and extending in to 


extinctions Ay >10, deep within the shielded inner lay¬ 
ers of a molecular cloud. For the comparisons below, we 
use the same chemical abundances, grain properties and 
rate coefficients as in Hollenbach et al. (2009) with the 
exception of the binding energies of O and OH on grain 
surfaces. For the binding energy of O, we adopt a char¬ 
acteristic value from this work, Ab(0) = 1800 K; for the 
binding energy of OH, we adopt the highest value from 
He & Vidali (2014), £’b(OH) = 4800 K, though the pre¬ 
cise value is not important so long as Ub(OH) > i?b( 0 ), 
as discussed below. Here we assume that F^des = Ab, 
i.e., that the binding (adsorption) and desorption ener¬ 
gies are the same, which is true if the processes are not 
activated. 

In the models, the abundances of species on dust grains 
are controlled by competition between adsorption onto 
grain surfaces, chemical processes which occur on the 
grain surfaces, and desorption processes—including pho¬ 
todesorption by FUV photons, thermal desorption from 
warm grains, and cosmic ray desorption. The chemistry 
of oxygen-bearing molecules was shown to be driven by 
grain-surface reactions. The chemistry of water, in par¬ 
ticular, begins with the adsorption of atomic oxygen onto 
a grain surface. If the grain encounters a hydrogen atom 
before the oxygen atom is desorbed from the surface, 
then OH ice is formed. In the model, it is assumed that 
H atoms are highly mobile on the grain surface can eas¬ 
ily find O in ice. If the OH is at least as tightly bound 
to the surface as the O, then OH will also abstract an 
H, forming water ice on the grain. This grain-surface ice 
is then the source, via photodissociation, of gas phase 
water near the cloud surface. The main result of the in¬ 
creased O binding energy described in this work is that 
O will stick to grains that are warmer, so that the grains 
can catalyze water formation in higher FUV fields and 
closer to the cloud surface. 

Following the discussion in §2.5 of Hollenbach et al. 
(2009), we find that the new value of the O binding 
energy should allow formation of water on grains for 
grain temperatures as high as ^ 50 K. In the models, 
the grain temperature in the unshielded cloud surface is 
Tgr ^ 12 G°'^ K. Thus grain surface formation of water 
should proceed efficiently near the cloud surface for FUV 
fields as high as Gq ^ 10 ^, and may also be efficient in 
regions with higher FUV fields deeper into the clouds 
where the grain temperature has dropped to ~ 50 K. 

To illustrate the effects of the increased binding ener¬ 
gies, we have computed models with Ub(0) = 800 K and 
with £'b(0) = 1800 K for two different combinations of 
n and Gq. In the first case, we compute the results for 
models with n = 10 ^cm“^ and Gq = 10 ^, conditions ap¬ 
propriate for a molecular cloud interface near the outer 
edge of an H II region. In Figure 5, we show the results 
for i?b(0) = 800 K and in Figure 6 , we show results 
for Ab(0) = 1800 K. The effects may be seen in the 
outer edge of the cloud, where the grain temperature is 
above 30 K. In the case of the old binding energy, O 
does not stick to grains and begin to form water ice un¬ 
til the grain temperature has dropped below ^25 K, at 
Ay ^ 0.5. For the higher O binding energy, O can stick 
to grains all the way out to the cloud surface, where the 
grain temperature is ^ 30 K. However, the strong FUV 
field is able to photodesorb water ice from the grains at 
a rate high enough to prevent a monolayer of ice from 
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Fig. 3.— TPD traces of mass 16 amu (blue) and mass 48 amu (green) after deposition of 0, 240, and 480 s of O/O 2 on 0 . 2 , 0.4, 0.6, and 
0.8 ML of O 3 pre-coated amorphous silicate. The heating ramp is 0.5 K s“^. The top row of each panel has the original TPD traces while 
in the bottom row the contribution of Oj fragmentation to the signal of mass 16 amu has been subtracted. 
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Fig. 4.— Fitting of the mass 16 amu TPD traces after subtraction 
in the top-right panel (4 min O 3 -|- 0, 4, 8 min O) of Figure 3 using 
= 1850 K. 


forming. Thus, although the water ice and gas phase 
water abundances in the outer regions of the cloud are 
several orders of magnitude higher than if the O binding 
energy were lower, the effect on the total column den¬ 
sities of water and other 0-bearing species is negligible. 
This model simply illustrates the effect of the binding 


energy on the ability of O to reside on grains. 


n=10* cm-s, G„=10=, E,(O) = BO0 K, E,(OH) = 4800 K 



Depth into the Cloud in Visual Magnitudes, 

Fig. 5.— Photodissociation region model for n = 10“^ cm~^ 
and Go = 10^. with O binding energy of 800 K. Left axis shows 
the abundances of 0 -bearing species as a function of cloud depth, 
while the right axis indicates the grain temperature. The dashed 
line shows the abundance of water ice on grains. 
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n=10* cm-=, G„=10', Ej(0)=lBOO K, Ej(0H)=4B00 K n=10» cm-=, G„=10‘, E,(0)=800 K, Ej(OH)=4aOO K 



Fig. 6 .— Same as Fig. 5 but with an O binding energy of 1800 K. 

A more extreme case can be seen in Figs. 7 and 8. Here 
we make the same comparison as in Figs. 5 and 6 but 
with n = 10® cm“^ and Gq = 10^, about the most ex¬ 
treme values expected in a dense star forming region like 
Orion or a dense nuclear starburst. Here there are two 
possible routes for water formation. In the outer layers 
of the cloud, the gas temperature can be several hun¬ 
dred degrees, high enough to drive the neutral-neutral 
reactions O -I- H 2 —> OH -|- H and OH -|- H 2 —H 2 O 
+ H, forming water directly in the gas phase. In addi¬ 
tion, O 2 can form from O -I- OH. These reactions are 
responsible for the broad water peak in the two figures 
at A\r ^ 0.3. Beyond this depth, the gas temperature 
is too low to drive these reactions, and water formation 
must depend on the grain surface mechanism. In the 
model with low binding energy, O stays off the grains 
beyond Ay ~1 and is driven into O 2 which reaches an 
abundance ~ 10“®. The fact that this was not observed 
is precisely what led Hollenbach et al. (2009) to suggest a 
higher binding energy. The effects of the higher binding 
energy on H 2 O and O 2 abundances can be seen in Fig. 
8 . Once the grain temperature drops below 60 K, the 
model with the high value of the O binding energy can 
begin to form water ice on grain surfaces. The result is 
an increase in the column density over the lower binding 
energy case of a factor 4, with a corresponding increase 
the in the intensity of optically thick but effectively thin 
water emission lines. The effect on the column of molec¬ 
ular oxygen, however, is in the other sense; with a higher 
value of the O binding energy, O is locked up in water 
ice and can not form O 2 by the O + OH reaction as it 
would if O were still in the gas phase. The result is a 
lower gas-phase column of O 2 than had been predicted 
for high Go cases in Hollenbach et al. (2009). In another 
paper, we will more extensively explore the effects of the 
new binding energies. 

4. SUMMARY 

Oxygen is the third most abundant element, and sig¬ 
nificant observational and theoretical work has been ded¬ 
icated to understanding its distribution in the interstellar 
medium. In this work, we explored the possibility that 
oxygen might reside on grains over a wider temperature 
range than previously accepted. In the first measurement 



Fig. 7.— Photodissociation region model for n = 10® cm~® 
and Go = 10^ with O binding energy of 800 K. Left axis shows 
the abundances of O-bearing species as a function of cloud depth, 
while the right axis indicates the grain temperature. The dashed 
line shows the abundance of water ice on grains. 


n=10* cm-s, G„=10‘, E,(0)=1800 K, E„(OH)=4800 K 



Fig. 8 .— Same as Fig. 7 but with O binding energy of 1800 K. 

of oxygen binding to dust grain analogs, we find that the 
values of the binding energy (1660 and 1850 K for O 
on, respectively, porous amorphous water ice and a bare 
amorphous silicate) are more than twice the previous es¬ 
timate (800 K, Tielens & Hagen 1982). Our values are 
near the ones that were suggested in order to reconcile 
observations with abundance constraints (1200 K and 
1600 K proposed by, respectively, Hollenbach et al. 2009; 
Melnick et al. 2012). Since the estimated value (800 K) 
has been widely used in simulations of the chemical evo¬ 
lution of ISM environments, we performed simulations 
with the old and new values using as a test the same 
ISM environment as in Hollenbach et al. (2009) and Mel¬ 
nick et al. (2012), i.e., a molecular cloud edge exposed to 
FUV. For a moderate illumination (Go=10^), the abun¬ 
dances of H 20 ice and H 20 gas increase by orders of mag¬ 
nitude, although it leaves the column density marginally 
affected. More interesting is the case of higher FUV il¬ 
lumination (Go = 10^), a case expected in a dense star 
forming region as in Orion. The higher gas temperature 
allows neutral-neutral reactions to form OH and H 2 O. At 
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Ay 1, the gas temperature has dropped enough that 
water ice formation occurs aplenty on surfaces of grains. 
While the column density of water increases, the column 
density of gaseous O 2 decreases because more oxygen is 
locked up in grains. 
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